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Description 

Computer-aided method for par allel calculation of the operating 
point of electrical circuits 



The computer-aided simulation of electrical circuits has attained 
increasing importance in the development of very large circuits, 
that is to say circuits having a very large number of elements. It 
is particularly in the development of computer chips having a 
10 multiplicity, for example several hundred thousand transistors, 
that serial processing for the determination of the circuit 
quantities by a computer has been found to be unusable because of 
the excessive time consumption. 

15 In WO 98/24039, therefore, it is proposed to partition a large 
circuit and to have the partitions processed by different 
computers in each case. 

In the calculation, the operating point, that is to say the 
20 potentials of all nodes is usually determined first as basis for 
further analyses such as, for example, transient or alternating- 
current analyses. 

For the parallel calculation, an inclement at ion of the Newton 
25 method is proposed in U. Wever, Q. Zheng et al.: "Domain 
Decomposition Methods for Circuit Simulation'' (Proceedings of the 
8th Workshop on Parallel and Distributed Simulation, PADS '94 
Edinburgh, Scotland, UK, pp. 183-186, July 1994) and in U. Wever, 
Q. Zheng: "Parallel Transient Analysis for Circuit Simulation" 
30 (Proceedings of the 29th Annual Hawaii International Conference on 
System Sciences, pp. 442-447, 1996) . The disadvantage is that 
convergence can only be achieved here when sufficiently good 
estimates of the operating point are available, due to poor 
convergence characteristics. As a rule, however, such good 
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estimates can be achieved with difficulty or not at all in the 
case of large circuits. 
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The present invention is, therefore, based on the object of 
creating a computer-aided method for parallel calculation of the 
operating point of electrical circuits which ensures simple, 
reliable and fast calculation pf the operating point of the 
5 circuit. 



10 



According to the invention, this object is achieved by the method 
as claimed in claim 1. 

In the method, the charging method, which has very good 
convergence characteristics, is- used for the parallel calculation 
of the operating point. 

j=i As described, for example, in H. Spiro: "Simulation of integrated 

15 circuits" (2nd edition, R. Oldenbourg Verlag Munich, Vienna 1990) , 

a 

the charging method uses the dynamic elements such as capacitances 
\Q and/or inductances existing in an electrical circuit for 

'ss? calculating the operating point via a pseudo-transient analysis. 

In this process, the dynamic elements are charged up step-by- step 
•i 20 to a value which, in principle, can have any magnitude. 



particularly to the value '^l", starting from a value zero. 



''~J In an embodiment according to the invention, dynamic elements such 

[zi as, for example, capacitances or inductances having predetermined 

^ 25 known values can be additionally installed at suitable places, for 

exan^le at or between a number of nodes. This makes it possible to 
avoid problems in the calculation due to the number of dynamic 
elements being too small. 

30 For this purpose, the dynamic elements additionally inserted into 
the circuit can be reduced step-by- step to zero from an initial 
value, for example tending to infinity, for which the operating 
point calculation is trivial, so that the original circuit is 
again simulated. 
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In a preferred embodiment of the invention, at least one node, 
that is to say a junction of at least two current paths, of the 
circuit is connected by means of a capacitance to a predetermined 
value having a predetermined potential. 

In a further development of the invention, however, a capacitance, 
the second terminal of which is in each case connected to a 
predetermined potential, for example to ground, can also be 
connected to each node via all partitions. This procedure has the 
advantage that the calculation of the operating point, i.e. of the 
respective potentials of the individual nodes, is trivial for the- 
circuit for an initial value for the capacitances which tend 
toward infinity, due to the equation, which is then explicit, 
being solved. Changing the value for the capacitances step-by- step 
then makes it possible to change the circuit simulation by 
suitable new selection of the value for the capacitances until a 
calculation of the operating point of the circuit is obtained for 
a value of the capacitances tending toward zero or almost toward 
zero . 

By also implementing this charging method for parallel calculation 
of an electrical circuit, a very large circuit having a 
multiplicity of transistors can be advantageously calculated in a 
fast and simple manner even with a smaller number of dynamic 
elements. Various procedures are conceivable for the respective 
new determination of the value for the capacitance with the aim of 
allowing this to go toward zero, the decision criterion being the 
degree of difficulty of calculating the operating point of the 
respective preceding step. 

The operating point obtained by solving a nonlinear equation for a 
particular value of the capacitances in each case can be solved, 
for example, iteratively by means of the Newton method. For 
choosing the next value for , 
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the number of iteration steps necessary for the solution can then 
be used until the value "drops below a predetermined value for "C" . 

Further advantageous embodiments of the invention are obtained 
from the subclaims. 

In the text which follows, the invention will be explained by 
means of an exen^lary embodiment shown in the drawing. In the 
drawing, the single figure shows a flowchart which represents the 
individual method steps of the method. 

According to the drawing, a simulation of an electrical circuit is 
the starting point 1, for example in the circuit description 
language SPICE. This is partitioned in a first processing process 
2, as described, for example, in WO 98/24039 as "Clustering 
method", so that individual partitions or parts of the circuits 
are obtained which can be calculated with the same degree of 
difficulty, if possible. 

Naturally, other applicable variants of partitioning are 
conceiveUble such as, for example, the ""ratio cut method" described 
in N. Frdhlich, B. Riess, U. Wever, Q. Zheng: *A new approach for 
parallel simulation of VLSI circuits on a transistor level" (IEEE 
Transaction on Circuits and Systems - I: Fundamental Theory and 
Applications, Vol. 45, No. 6, June 1998, pages 601 to 613), or 
even an arbitrary "manual" division. 

In a further process step 3, a grounded capacitance is added to 
each node of the circuit, that is to say to a junction of at least 
two conductors or current paths, respectively. Naturally, it is 
also conceivable to add a capacitance, the second terminal of 
which is connected to a predetermined potential, in which case 
both potentials and values for the capacitances of each node 
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can differ. For reasons of efficiency, stsibility and economy of 
calculation, each capacitance is connected to the same potential 
with its second terminal, for example ground, an identically high 
value CO also being selected for all capacitances in step 4 for 
the same reasons. 

For this value CO, the operating point is then calculated for each 
partition or part -circuit in a further step 6, in which the 
required coupling values, that is to say the values for the 
coupling points or interfaces of adjacent partitions being 
exchanged and included in the calculation of the operating points 
of adjacent partitions. 

In this process, a partition which will be called "master" in the 
text which follows can advantageously take over the control of the 
charging process for reasons of efficiency. The master then 
determines the initial value CO for the capacitances, and it is 
also conceivable to predetermine the initial value externally, for 
example by the user. This value CO is then transferred to all 
other partitions, called slaves in the text which follows. 
Following this, the operating point is calculated both in the 
master partition and in all slave partitions, during which 
process, naturally, the full source vector, which represents the 
energy sources existing in the circuit, is present at the circuit. 

Starting from an overall problem to be solved 



if • • ' f 




where m is the unknowns of the master, 
p is the number of partitions, 
Si is the unknowns of partition i. 
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the charging process leads to the system of differential 
equations. 
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f(x,t) + D^ = 0 

where 

D = diag (C, ... C, 0, 0) and 

t is the time. 

This system of differential equations can be solved, for example, 
by means of the implicit Euler method which leads to the. non- 
linear equation 

where is the solution at time t*" for k = 0, 1, etc and h is the 
step t^*^ - t*'. This non-linear equation can be solved, for example, 
iteratively by means of the Newton method 

for n = 0, 1, etc. 

In the parallel calculation, the master then only calculates the 
unknowns ''m", the unknowns "si" representing fixed areas for the 
master. To be able to calculate a solution of this, the slaves or 

k+l 

slave partitions calculate the corrected values Sj in each 
iteration step of the master and report these to the master. For 
this purpose, the slaves must solve the system of non- linear 
equations. 
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[sr.mr.r)+^D,((sr'.mr')-(sNm^))=o 



where fi and Di are the corresponding equations and matrices for 
the slave '^i* . This non-linear equation can also be calculated by 
means of the Newton method and not all iteration steps have to be 
necessarily performed to convergence. 

The results of the slaves are then inserted into the system 
equations of the master whereupon the master can calculate the 
expression tn^^j^ . 

In a next step 7, a new value ''Cnew" for **C" is determined' by the 
master and the slaves can also make suggestions. To achieve the 
aim, namely a value for "C" tending to zero or, respectively, , less 
than a predetermined value e, and thus a calculation of the 
original circuit, the choice of the new value '^Cnew" tor «C" is 
made dependent on the difficulty of calculating the preceding step 
with the value '^Coid" for . 

In this context, various procedures are conceivable, for example 
the choice of ^Cnew" in dependence on an analysis of the number of 
iteration steps which the master needed for calculating the non- 
linear system of equations for the preceding value ^^Coid" for '^C : 




2 IS n < 

isHj <n<n2 
•2 isnj >n 
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where ''n" is the number of iterations of the master and *'ni" "n2" 
are the parameters predetermined by the user. Naturally, instead 
of halving or doubling ''Coid" / respectively, other strategies for 
reducing or enlarging "Coid" are also possible. Additionally, the 
number of iteration steps needed by the slaves or slave partitions 
for solving their non-linear system of equations can also be taken 
into consideration, for example through the choice of 

n = max n,^ni,slavel,...,^npSlavep 

V i=l i=l J 

where ni, siavej is the number of iteration steps of slave j, during 
the i-th iteration of the master. It must be noted in this context 
that the choice of a value for only affects the efficiency of 
the method not the operating point itself. 

After a value for *vc" of less than or equal to a predetermined 
value e has been reached, the calculation is terminated at branch 
5 and the value for C can be set to **.0" in a last step. The 
operating point of the original circuit is thus found as result 8 
and can be output via output units such as for example, a screen, 
printer or the like and/or stored in a memory as basis for further 
analyses of the circuit. 

Due to the method according to the invention, an operating point 
of a very large electrical circuit can be advantageously 
calculated in parallel by a multiplicity of computers or 
processors and the disadvantage of known parallel types of 
calculation, namely lack of convergence with unfavorable starting 
values can be avoided. Due to the parallel calculation, the 
expenditure for calculation with regard to iteration steps and 
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setting a value for "C" step by step, can be kept within tolerable 
limits due to the distribution to a number of processors. 

The operating point thus calculated is then used as basis for 
further analyses, for example the alternating- current analysis of 
a circuit. 



